Quantum impurity models of noisy qubits 
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We demonstrate that the problem of coupled two-level systems ("qubits") which are also subject 
to a generic (sub)Ohmic dissipative environment belongs to the same class of models as those 
describing (non)magnetic impurities embedded in strongly correlated systems. A further insight into 
the generalized single- and two-impurity Bose/Fermi Kondo models enables one to make specific 
recommendations towards a systematic engineering of highly coherent multi-qubit assemblies for 
potential applications in quantum information processing. 



The studies of open quantum systems have long been 
focused on the behavior of a single two-level system (an 
effective spin-1/2 or, as it is nowadays often referred to, 
a qubit) in a dissipative environment, for this hallmark 
problem appears to describe a variety of physically dif- 
ferent (yet, formally related) situations [1]. 

More recently, the attention has been drawn to this 
problem's further extension where N coupled two-level 
systems would be subject to a dissipative (possibly, non- 
uniform) environment. This interest was largely moti- 
vated by the relevance of the problem in question for 
practical implementations of the ideas of quantum infor- 
mation processing. 

However, despite the previous (mostly numerical and 
largely limited to weak coupling) analyses of the problem 
of iV = 2 noisy qubits [2], not much of a consideration 
has been given to the assessment of a possibility of reduc- 
ing the environmentally induced decoherence by means 
of properly choosing the parameters of individual qubits 
and/or by virtue of permanent pairwise qubit interac- 
tions which (whether desired or not) would be unavoid- 
ably present in any realistic setup. 

As one step towards developing a systematic (as op- 
posed to heuristic) approach to addressing this kind of 
issues, in the present work we explore a formal connection 
between the spin-boson model of interacting qubits cou- 
pled to a generic (multi-component and/or (sub)Ohmic) 
bosonic bath and some of the recently studied models of 
(non)magnetic impurities in strongly correlated electron 
systems. 

Such analysis would be incomplete without investigat- 
ing the intermediate-to-strong coupling regime where the 
customary Fermi's Golden rule-based estimates fail, and 
a more elaborate approach is needed. While, at the first 
sight, a strong-coupling behavior may not seem to be 
immediately relevant to the previously proposed designs 
of a practical quantum register, it should be noted that 
experimentally measured decoherence rates routinely ex- 
ceed their best estimates obtained in the framework of 
the weak-coupling Bloch-Redfield and related approxi- 
mations. Moreover, even in the case of a single qubit, the 
recent non-perturbative analyses of the strong coupling 
behavior as well as the effects of non-Markovian and/or 
structured environments have already resulted in a num- 



ber of potentially interesting implications for quantum 
computing [3]. 

Regardless of the qubits' physical make-up, the dy- 
namics of N such species is usually described by the spin- 
boson Hamiltonian 
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The random field /if = J2k ^''''''KiK + b-k)/y^ repre- 
sents a generic multi-component bosonic bath composed 
of Z?-dimensional propagating modes with the dispersion 
LUk = vk and described by the correlation function (here- 
after i? = |j — j| is a distance between the and j 
qubits and A is an upper cutoff of order the bath's band- 
width) 
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The variable parameter < e < 1 controls the bath's 
spectral density p{uj) = X]/c(l-^fel^/^'=)'^(^ ~ ^k) oc 
A'^u!-^~'^, thus allowing one to study the entire range of 
different (sub)Ohmic environments within the same uni- 
fying framework. 

Recently, a new reincarnation of the Hamiltonian (1) 
has emerged under the name of the Bose Kondo model in 
the theory of (non)magnetic impurities in strongly corre- 
lated systems, such as heavy fermions, Mott insulators, 
and high- Tc cuprates [4]. 

In the case of a single qubit {N = 1), the Hamilton (1) 
gives rise to the renormalization group (RG) equations 
which, to the first order in the dimensional regulariza- 
tion parameter e — 3 — D <^ 1, read 
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In these equations which gcncrahzc the analysis of Ref. 
[4] to the case of a finite uniform field B ^ Q, I = 
InA/w is the energy-dependent renormalization scale. 
The longitudinal g\\ = {■yB)'^/B'^ and the transverse g±^a 
iJ2a=i 2 .9-L,a = 7^ ^ {jB)^ / B"^) qubit-bath couplings are 
defined with respect to the direction oi B = (0, 0, A). 

Eqs.(3) feature a variety of fixed points that can be 
classified according to the number of effectively equal 
components of the vector g. Namely, we find unstable 
SU{2)- (5_L,a = 511) and XY- {gj_^a. = 9, 9\\ = 0) invariant 
fixed points, while the only stable one is a pair of Ising 
fixed points (gy = g±^i = or (jy = g± 2 = 0). Thus, in 
a striking contrast with the conventional (Fermi) Kondo 
model [1], an anisotropy of the bosonic couplings does 
not renormalize away, but instead tends to increase. 

In the Ising case, the coupling to a one-component 
Ohmic bath gives rise to the conventional Kostcrlitz- 
Thouless (KT) transition that occurs at 57 = 1, A/ = 
and separates the regime where a sufficiently small ini- 
tial A continues to decrease (g > f) from that where 
it instead grows {g < 1). In the language of the Fermi 
Kondo model (see below), these two regimes correspond 
to the ferromagnetic (FM) Kondo effect and the usual 
antiferromagnetic (AFM) Kondo screening, respectively. 

The behavior of the running RG variable A(Z) should 
not, of course, be confused with that of the renormal- 
ized splitting between the qubit's energy levels A{1) = 
A(i)e~' which still decreases (albeit, in different ways) 
in both cases. In the extensively studied AFM regime of 
the Ohmic problem, the latter is given by the well-known 
solution 

A, ~ A(Ao/A)i/i-» (4) 

of the self-consistent equation A, = A(/ = In A/ A*) [1]. 

In the case of a sub-Ohmic (e > 0) bath and in the 
vicinity of the Ising fixed point, the RG flow decribc^d by 
Eqs.(3) is characterized by a constantly increasing effec- 
tive coupling g{l) and a vanishing (possibly, after some 
temporary increase) A(Z), provided that its bare value 
Aq falls below the separatrix which connects the unsta- 
ble fixed point at = 1, A/ = e^/^ to the trivial one at 
the origin of the g — A plane. 

Such a behavior signals a total loss of coherence and 
a complete qubit's localization in one of the two de- 
generate states (overdamped regime), in agreement with 
the customary expectations [1]. However, for Aq ^ 
A[2ge^^'^~^]^^'^ the solution of Eqs.(3) shows that coher- 
ent damped oscillations with a frequency A* ~ A(2gy^'^ 
might still be possible, thus supporting the earlier pre- 
diction made in Ref. [5]. 

In contrast to the Ising case, the coupling to a multi- 
component ("non-abelian") Ohmic bath may give rise to 
a markedly different behavior. In both cases of the XY- 
and SU (2)-symmetrical couplings, our solution of Eqs.(3) 
shows that the standard KT fixed point is absent, and the 



renormalized level splitting continues to grow (although 
it may temporarily descrease at g{l) ~ 1), thus enabling 
the qubit to avoid localization. 

Furthermore, in the "weakly sub-Ohmic" (0 < e <C f ) 
case a new unstable fixed point emerges at gxY = e 
{gsu(2) = e/2) for the XY- {SU{2)-) symmetrical cou- 
plings and A.xY,su{2) = 0. Notably, this non-trivial fixed 
point governs the regime of small g{l) and which is 
of a primary importance for quantum computing-related 
applications of Eq.(l). 

In the latter case, contrary to the naive expectation 
drawn from the weak-coupling analysis, the coupling g{l) 
may first increase, before succumbing to the continuing 
growth of A (I) and reverting to the opposite, decreas- 
ing, behavior, provided that the initial values of the RG 
variables satisfy the condition max[gQ, Aq] ^ e. 

Besides, Eqs.(3) indicate that for a strongly sub- 
Ohmic bath (e > e, = 1/2) yet another unstable XY- 
symmetrical fixed point can emerge at gxY,ii = 1/2 and 
AxY.ii = (e ^ e*)^^^, while the fixed point at gxYj = e 
and AxY,i = would then become stable. Considering 
the fact that the corresponding threshold value e* is quite 
large, this prediction may turn out to be spurious. Nev- 
ertheless, from a general viewpoint, it is conceivable that 
even a multi-component (yet, sufficiently strongly sub- 
Ohmic) bath could still be capable of causing a complete 
localization of the qubit for small Aq. 

The above predictions for the RG flow of the couplings 
and the effective level splitting have a direct bearing on 
the properties of such observables as the qubits' correla- 
tion functions < {S^{t), Sj{0)} > and the corresponding 
dynamical spectral functions 

<n\S-\m><m\S-\n> 

Xij (w) = Im2_^{pnn - Pmm)— p , p , -p 

(5) 

given by the sums over renormalized and broadened en- 
ergy levels |n > of the interacting noiseless qubits which 
are weighted with the diagonal elements of the equilib- 
rium density matrix /9„„ = \n >< n\. 

The analysis of Eq.(5) reveals that for non- interacting 
qubits the coherent (single-qubit) behavior manifests it- 
self through the presence of inelastic peaks at w = ±A*/2 
in the normalized spectral density Xij{^)/^^~'^ cx: Sij 
[3,5]. The width of these peaks is controlled by the de- 
coherence rate 

r (X 5*p(A,) coth(A*/2T) (6) 

which increases as F ~ g^T{K/ A^Y with increasing 
and decreasing A* < T. As the qubits lose their coher- 
ence, the spectral weight gets transferred from the inelas- 
tic peaks to the small energies (w ~ 0), which behavior 
constitutes the decohering effect of the environment. The 
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onset of complete localization is usually preceeded by ex- 
ponential relaxation which separates the former regime 
from that of the partly coherent damped oscillations [1]. 

It is worth reminding that different probes may sig- 
nal the loss of coherence at different critical couplings. 
In the extensively studied Ohmic case, such a difference 
is cxamplificd by the juxtaposition of the singlc-qubit 
average < S^{t) > and the auto-correlation function 
Xii{w) cothu}/2T whose coherent peaks get washed out 
at, respectively, g = 1/2 (Toulouse point) [1] and g = 1/3 
[6], both critical values being lower than the beforemen- 
tioned estimate gi = 1. 

In light of the above mentioned possibility of a non- 
monotonic behavior of g{l) near the SU{2)- and XY- 
symmetrical fixed points at small go and Aq, the initial 
increase of the coupling strength may give rise to the sit- 
uation where the inelastic peaks in the spectral function 
(5) might temporarily become rather broad before they 
can narrow down at still lower energies. 

The above findings strongly suggest that by operat- 
ing a noisy qubit in the vicinity of one of the XY- and 
S't/(2)-symmetrical fixed points one can better retain its 
quantum coherence. 

One could argue, however, that for non-interacting 
qubits the intrinsic instability of the model (1) towards 
developing the Ising-like anisotropy may hinder any pos- 
sibility of taking advantage of the greater robustness of 
quantum coherence found at the high-symmetry fixed 
points. To this end, below we show that the loss of co- 
herence caused by the outward flow of the single-qubit 
RG trajectory away from a desired operating point can 
be thwarted by inter-qubit couplings which provide for 
an extra protection against decoherence. 

In fact, some exchange-like interaction between the 
qubits is necessarily generated by the qubit-bath cou- 
plings themselves. The instantaneous part of this (gen- 
erally, retarded and FM-likc) interaction which arises in 
the course of integrating over the bosonic modes down to 
the energy scale u ~ v/R 
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combines together with the direct inter-qubit coupling 
into the effective parameters introduced in Eq(l). 

It can be easily seen that a sufficiently strong FM ex- 
change (/ — oo) forces a pair of qubits into a triplet 
state which then follows an effective S = 1 Bose Kondo 
scenario. In the opposite, strongly AFM, limit (J +Qo) 
the two qubits get locked into a singlet state, which pre- 
vents them from any unwanted entanglement with the 
environment. According to Eq.(7), however, the latter 
regime can not be attained in the absence of a sufficiently 
strongly AFM direct coupling between the qubits. 

At intermediate values, / provides a cutoff for the RG 
flow which now terminates at = In A/|J| before reach- 



ing the strong coupling limit, thereby resulting in a larger 
effective A* = A(//) that determines the position and the 
width of the coherent peaks. 

In the case of the 3D Ohmic bath, a further insight can 
be gained from the previously established correspondence 
between the N = 2 Hamitonian (1) and the anisotropic 
two-impurity Fermi Kondo (TIFK) model 
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where the auto-correlation function of the ID spin- 
1/2 fermions with the Fermi momentum kp is < 
c|(i)cj(0) >(x 6ij/t, while the Kondo couphngs 
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between S'^ = ± S"! and the even, odd, and mixed 
bilinear combinations of the fermion fields Ci,2 = (ce ± 
Co)/V2 at the two qubits' locations are given in terms of 
the exchange constants of the single-impurity anisotropic 
Fermi Kondo model [7] (which, in turn, are related to 
the parameters of the N = 1, e = Hamiltonian (1): 
5 = (1 - (2/7r) tan-i(7rjll /4A))2 and A = [1]). 

In contrast to its bosonic counterpart (7), the RKKY- 
like inter-qubit coupling mediated by the fermionic 
bath behaves as Ip cx {,P /A){{2kFR) cos{2kpR) — 
sm.{2kFR))/{2kFR)^- Therefore, unless one chooses to 
neglect any environmentally induced contributions to the 
exchange coupling / altogether, the strict correspondence 
between TIKM and the = 1 model (1) can only hold 
for small inter-qubit separations, R <C l/kp (see below). 

In the presence of the particle-hole symmetry, the two- 
channel TIKM (8) possesses an unstable fixed point at a 
critical (AFM) value 7* which is set by the single-qubit 
Kondo scale T^~^^'^\ This critical point, which in the 
S'?7(2)-symmetrical case (r^^=i/2) = Aexp(-A/J)) oc- 
curs at /* « 2.5T^~^^^\ gets replaced by a crossover, if 
the particle- hole symmetry is broken [7]. 

In the FM regime (/ < -T^j^-''^'), the Kondo screen- 
ing of the composite spin = 1 is characterized by 
the scale Tk{I) ~ r^s=i/2)(2.^s=i/2)^|^|^^ ^j^.^j^ 

polates between T^^=i/2) and T|=i = Tk{I < -A). 

The non-universal exponent r] depends on the exchange 
anisotropy and assumes the value 77 = 2 in the SU{2)- 

Aexp(-2A/J)). 



invariant case (T^ 
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The above behavior pertains to the two-channel model 
which describes a pair of well separated qubits {R ^ 
l/kp), whereas at small separations the term in (8) pro- 
portional to J° <C J+ ~ can be disregarded, and 
Eq.(8) reduces to the single-channel TIFK model formu- 
lated solely in terms of the even combination Cg of the 
fermion orbitals which is coupled to the total spin 5+. 
Remarkably, the two- and single-channel TIFK models 
appear to describe the two opposite limits which corre- 
spond to independent and collective decoherence, respec- 
tively. 

For a pair of identical qubits (Ji = J2), the single- 
channel TIKM undergoes a simple first order transition 
(level crossing), as / is tuned past its critical value /*, 
from the underscreened S = 1 model at / < /* to the two- 
qubit singlet state which is completely decoupled from 
the bath at / > [8]. 

If, however, the qubits are different (Ji 7^ J2), then, 
as / increases past the more strongly coupled qubit 
first gets screened at T^^ = Aexp(— A/maa;[Ji, J2]), fol- 
lowed by the KT-type transition (which now occurs re- 
gardless of the presence of the particle-hole symmetry) 
at T^"' max[I - h.T^^^ cxp{-TP/I - Q]. 

The above conclusions, which were drawn in the case of 
the 5J7(2)-invariant couplings, apply not only to qubits 
represented by physical (electron or nuclear) spins but 
also to the effective qubit operators which represent, e.g., 
electron states in two-level single or lateral/ vertical dou- 
ble quantum dots. Such operators have been routinely 
introduced in the Anderson or resonant level models 
where the genuinely S'?7(2)-invariant exchange couplings 
are generated from the electron tunneling amplitudes by 
the Schrieffer- Wolff transformation. 

A comparison between several different kinds of the 
Ising-type qubit interactions studied in the context of 
the N = 2 problem [2] indicates that a pair of qubits can 
better retain their coherence if the interaction term com- 
mutes with the qubit-bath coupling operator. In this 
regard, the Heisenberg exchange, which not only com- 
mutes with the overall qubit-bath coupling in the case of 
collective decoherence but also provides for the strongest 
lifting of the singlet-triplet degeneracy, might offer the 
best available option for reducing decoherence. 

We note, in passing, that spin-rotationally invari- 
ant couplings may not be easily realizable for some of 
the pseudo-spin qubits. Among such examples are the 
Josephson junction qubits whose Hamiltonians, while be- 
ing potentially well-controllable, generally would have no 
particular symmetries at all. Therefore, it is conceiv- 
able that future designs of practical solid-state qubits will 
need to optimize between the somewhat contradictory re- 
quirements of efficient control and robust coherence. 

The formal analogy between the problems of noisy 
qubits and quantum impurities established in this Letter 
can be pursued even further. In particular, the B ^ 



counterpart of the previously studied Fermi-Bose Kondo 
model [4] can describe a system of qubits coupled to 
one Ohmic (represented by a fermionic) and the other 
sub-Ohmic (bosonic) bath. Such a situation occurs, e.g.. 
in the Josephson qubits where the effects of both the 
Nyquist (e = 0) and 1// (e = 1, if treated as approxi- 
mately Gaussian) noises which are caused, respectively, 
by fluctuating currents and background charges often 
need to be considered on equal footing. 

The fact that the only stable fixed points of the mixed 
Fermi-Bose model are the pure Fermi and pure Bose ones 
[4] suggests a possibility of effectively blocking off a more 
harmful type of coupling, while tackling the remaining 
one with, e.g., error-correction techniques. 

To summarize, in the present work we exploited a for- 
mal similarity between the problem of interacting qubits 
in (sub) Ohmic dissipative environments and the recently 
studied models of quantum impurities in strongly corre- 
lated systems. This insight enabled us to formulate a 
number of concrete recommendations for achieving the 
optimal (coherence-wise) regime for operating a noisy 
multi-qubit quantum register. 

Firstly, we investigated the RG properties of the 
(sub)Ohmic single-qubit model (1) in the vicinity of its 
fixed points of different symmetries. In the course of this 
analysis, we discovered that for non- interacting qubits 
and a "weakly-sub- Ohmic" bath (0 < e < e*), the best 
possible conditions for preserving coherence would re- 
quire the most spin rotationally-symmetrical qubit-bath 
couphngs 7". 

Secondly, we ascertained the possible benefits of per- 
manent interactions /" between the elements of a multi- 
qubit array for attaining the most coherence-friendly 
regime and protecting the qubits' quantum memory dur- 
ing the idling periods between consecutive gate opera- 
tions. Specifically, we found that it would be advanta- 
geous to time the (preferably, SU (2)-symmetrical) pair- 
wise intcr-qubit coupling / close to (yet, smaller than) 
the critical value /, of the TIFK model, thereby causing 
an incipient singlet formation and concomitant quench- 
ing of the Kondo screening. 

Lastly, the Kondo physics- conscious approach to the 
engineering of robust qubit Hamiltonians might prove in- 
strumental in solid state implementations of quantum in- 
formation processing where, unlike in liquid-state NMR 
or trapped-ion designs, such active noise suppression 
techniques as dynamical decoupling/recoupling schemes 
may not be readily available. Therefore, we believe that 
this work will further spur the ongoing cross-fertilization 
between the developing theory of quantum information 
and such well-established topics in materials theory as 
Kondo physics of heavy fermions and quantum dots. 

The author acknowledges valuable discussions with S. 
Girvin, A. MacDonald and S. Das Sarma. This research 
was supported by ARO under Contract DAAD19-02-1- 
0049 and NSF under Grant DMR-0071362. 
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Note added: In the case n = 1 and e > 0, Eq.(2) 
of this work (first made available as a preprint [cond- 
mat/0305510]) and the conclusions regarding the possi- 
bility of a delocalization transition for a sufficiently large 
A were recently confirmed in Ref. [9]. 
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